Anomalous Hall effect in superconductors with spin-orbit interaction 
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We calculate the anomalous Hall conductance of superconductors with spin-orbit interaction and 
with either uniform or local magnetization. In the first case we consider a uniform ferromagnetic 
ordering in a spin triplet superconductor, while in the second case we consider a conventional s-wave 
spin singlet superconductor with a magnetic impurity (or a diluted set of magnetic impurities). In 
the latter case we show that the anomalous Hall conductance can be used to track the quantum 
phase transition, that occurs when the spin coupling between the impurity and electronic spin 
density exceeds a certain critical value. In both cases we find that for large spin-orbit coupling the 
superconductivity is destroyed and the Hall conductance oscillates strongly. 

I. INTRODUCTION 

Anomalous Hall effect (AHE) was observed in metallic ferromagnets long time ago as a Hall current generated by 
electric field in the absence of external magnetic field . Since then several physical mechanisms of the AHE 
have been proposed, related to the side-jump and skew scattering from impurities PHS], inhomogeneous internal 
magnetization [110, internal spin-orbit interaction [S], and topology of electron energy bands jHllin]- Theory of AHE 
attracted much attention recently [3l [TTHT3| because it reveals some very unusual properties of solids such as the 
existence of monopoles in the momentum space or generation of topological gauge fields. 

One of the most intriguing models of AHE is the one based on an intrinsic mechanism [^110] related to the nontrivial 
topology of electron energy bands. In frame of this mechanism, the main contribution to AHE is due to electron states 
well below the Fermi energy [14'. The simplest model, in which this mechanism of AHE can be realized, is the model 
of a magnetized two-dimensional electron gas with Rashba spin-orbit (SO) interaction [15j. Unfortunately, it turns 
out that if the system is in the metallic state, i.e., if there is no gap at the Fermi surface, then the contribution of 
electron states at the Fermi surface can totally compensate the other contributions so that the resulting off-diagonal 
conductivity is zero |16| . In the opposite case, when the chemical potential lies in the gap, the anomalous Hall 
conductivity a^y is nonzero and quantized in units of /h. The theory of quantized AHE is quite similar to the 
theory of integer quantum Hall effect, where the gap is due to the Landau quantization in a strong magnetic field 

In this work we consider a 2D electron gas with nonzero magnetization and Rashba SO interaction. Such a model 
was used earlier for a description of intrinsic AHE |TB]. However, we calculate the AHE in the case when the electron 
system is additionally superconducting. The superconductivy produces a gap at the Fermi level, suppressing the 
contribution to AHE from the Fermi surface. Thus, one can expect that only the filled electronic states below the 
gap contribute to the AHE. The possibility of AHE in superconductors has been already considered in the case of 
ferromagnet-superconductor double tunnel junctions [19J, where side jump and/or skew scattering from impurities 
have been assumed as possible physical mechanisms responsible for the effect. This is, however, essentially different 
from our model, where we consider the intrinsic mechanism of AHE. Since in a superconductor charge is not conserved 
due to the particle-hole mixture, we do not expect any quantization of the anomalous Hall conductance. This was 
already shown for the usual Hall conductance in conventional superconductors in very high magnetic fields, where the 
Landau level description is appropriate |20| . 

Various materials are known to show the coexistence of ferromagnetism and superconductivity |21fl28| and, in 
particular, the presence of spin-orbit interaction due to the lack of spatial inversion symmetry |29H32| . We note, 
that the possibility of magnetoelectric effects in non-centrosymmetric superconductors was predicted already long 
time ago [53", where it was shown that a supercurrent should induce a spin polarization and reversely a Zeeman- 
like term should induce a supercurrent [34 as a result of strong spin-orbit interaction. Other effects due to the 
interplay of ferromagnetism and superconductivity have also been considered |35If57] . Recently, the interplay between 
superconductivity, magnetism and spin-orbit interaction (or topological insulators |38l I39| ) has received additional 
attention due to the possibility of Majorana edge states in a finite system or inside superconducting vortices [40 42J, 
with its possible applications in topological quantum computation. Moreover, the coexistence of magnetism and 
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superconductivity turned out to be interesting also from the point of view of possible applications in spintronics 

[13111]. 

In this work we consider in section II a spin triplet superconductor while in section III a conventional superconductor 
with a magnetic impurity |45j . In both cases we analyze the influence of Rashba spin-orbit interaction. In the first 
case the magnetization is due to a ferromagnetic order, whereas in the second case the system is locally polarized 
by a magnetic impurity. The latter situation may also be achieved when considering a superconducting film with 
a magnetic dot justaposed. It has been shown before that if the coupling between the magnetic impurity and the 
spin density of conduction electrons is strong enough, the system becomes magnetized through a first order quantum 
phase transition H7| that leads to discontinuities in various physical quantities |15] . In both cases we calculate the 
anomalous Hall conductance. We show that the Hall conductance of a superconductor with the magnetic impurity 
can be used to reveal the quantum phase transition. Finally, we conclude with section IV. 



II. AHE IN A TRIPLET SUPERCONDUCTOR 



We consider first a superconductor with a uniform magnetization. Since magnetism and superconductivity compete, 
a spin singlet superconductor is not stable due to Cooper pair breaking. Therefore, we consider a spin triplet super- 
conductor, where magnetism and superconductivity can coexist. The system is described by the tight-binding model 
in two dimensions, to which we add a superconducting pairing term with the appropriate symmetry. Additionally, 
we also include the Rashba spin-orbit term [TS] , which is generally allowed in non-centrosymmetric materials. Due to 
the spin-orbit term, a spin singlet component, Aj, is generally induced and therefore there is a pairing mixture in the 
system [49 . 

We write the electron operators, V'fc o-' terms of the Bogoliubov operators, ^, as 



(1) 



where fc, n label the eigenstates of the system. The wave functions and energy eigenvalues satisfy the Bogoliubov 
de Gennes equations [50) . which can be written as 
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— is the kinetic part, where t denotes the hopping parameter set in the following as 

the energy scale, t — 1, ep is the chemical potential, chosen in the following as ep = —1, fc is a wave vector in the xy 
plane, and we have taken the lattice constant to be unity, a = 1. Furthermore, in Eq.(2) is the magnetization, in 
energy units, along the z direction, while the vector d = {dx, dy,dz) is the vector representation of the superconducting 
pairing (p-wave). Finally, the Rashba spin-orbit term is written as Hfj ^ s ■ a ~ a {smkydx — ^vnkx<Jy), where a is 
measured in the energy units, and Ux, Uy are the Pauli matrices. 
The pairing matrix can be written as [ST] 
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Thus, we can write dx = [A^^^ — A^^) /2, dy ~ —i (A^^ + A^^^) /2, and dz = A^^, while the vector q ~ id x d* 



is given by Qx = Re (A;,; + A^,^) A* 
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vector vanishes, the pairing is called unitary. We have verified that considering the s-wave component has generally 
a very small effect on our results, and therefore we assume Ag = in the following. 
The energy eigenvalues of Eq. [2] can be written (for A^ = 0) as 
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and «!, 0^2 = ±. 

In the normal phase (c? — 0), the spin-orbit coupHng Hfts the spin degeneracy of the energy bands in the tight- 
binding model, except at fc = (0,0), (7r,7r) and (0,7r) (and equivalent points). These remaining degeneracies are lifted 
when including the magnetization. This is shown in Fig. [l] where the two energy bands are shown as a function of 
momentum for A^o = a/2 — 2 and various values of h^- As can be seen from Eq. |4] the lowest band is gapless at the 
points where 



{s-s + hl)+el = 2^{s-s + hl)el. (6) 

In a general case (d 7^ 0), the lowest band has gapless points that are solutions of the equation zi = 2^/z2, which 
yields 



d- d+ s- s + 4 + 2^d- s)' + (e| + dl)is- s + hi). (7) 



Thus, in the superconducting phase the system is generally gapped. In particular, without the spin-orbit interaction 
the gapless points are obtained by d • d -I- e| = which implies particular values for the chemical potential. 

The charge current along a link in the lattice can be obtained by adding a vector potential to the kinetic and 
spin-orbit terms and taking a functional derivative of the Hamiltonian with respect to the vector potential |50l I52| , 
or through its definition in the charge continuity equation [53J. The zero-momentum charge current in the = x,y 
direction can be written as 



— T 
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is a velocity matrix operator [M]. Here ?7^^ = cos (jt ■ and rj'i = sin (jt ■ (5^^ , where (5^ is a vector displacement 

(in units of the lattice constant) between nearest-neighbors along the fi direction. In turn, / is the 2x2 unit matrix. 

The Hall conductance can be now calculated using a Kubo like formula [55] . which in the limit of uniform and 
stationary current, if — > and cj — > 0, is given by 

/ ^ . ^ \ ^\ ^\ ^\ ^ fnk f m.k 



Q,;9 7,(5 n,m ^ ri.k m,k ' 
fe,/3^|7,5""(^' ")*«"(^' '^)«™(^' /^)"™(^' 7)* 

where TV is the number of sites and ^ is the Fermi function for the state described by n and k. In the normal 
phase the wave functions u and v are decoupled. The presence of superconducting pairing mixes the particle and hole 
character and, as already mentioned above, charge is no longer a good quantum number. The results for the Hall 
conductance depend then on the choice of the pairing matrix [511 156| . 

Let us now assume that the pairing amplitude is a free parameter. This describes the situations where supercon- 
ductivity is induced by proximity and therefore no self-consistent solution is implied. This also applies to a situation 
where a^y is measured on a normal sample in which superconductivity pairing exists due to proximity effect in the 
presence of a nearby triplet superconductor. We consider both unitary and nonunitary cases. Then we consider the 
case, where the pairing amplitude is determined by solving the Bogoliubov - de Gennes equations self-consistently. In 
the latter case we consider a nonunitary situation, for which the amplitudes A-^^^ and A^^^ are real, to simplify. This 
in turn implies that dy is imaginary. In all cases we take A^^ = (d^ = 0), which means that only the component 
may be nonvanishing. 

In Fig. [2] the anomalous Hall conductance in the normal phase (zero pairing amplitude) is plotted as a function of 
the magnetization hz and spin-orbit coupling A^o- The Hall conductance vanishes if either the magnetization or the 
spin-orbit coupling vanishes. Then, the absolute value of the Hall conductance increases as either parameter increases. 
Dependence on is more complex, as the conductance reaches a minimum around /12 = 1 ~ —ep, which shifts if we 
change the chemical potential. The minimum in the Hall conductance as a function of the magnetization hz (keeping 
the spin-orbit constant, for instance Xso = 2) is associated with the gaplessness of the spectrum at the point (0,7r) 
and the equivalent points (see also Fig. [ij. 

Now, we consider the superconducting phase. Since the spin-orbit coupling renders the type of pairing undefined 
(with the mixture of spin triplet and spin singlet pairings), the strength of the triplet pairing is expected to be 
weakened in comparison to the same superconductor with a vanishing spin-orbit coupling. However, it has been 
shown before [57] that the amplitude of the triplet pairing is not affected by the spin-orbit term when the vector d 
is parallel to the spin-orbit vector s. We have found that this pairing choice leads to results for the anomalous Hall 
conductance, that are very similar to those for the Hall conductance in the normal phase. This indicates that for this 
particular case, the superconducting order does not change significantly the Hall conductance, and therefore we do 
not show the corresponding results. We have also considered other choices of pairing, for which the vector d is not 
parallel to the spin-orbit vector s. We have considered both unitary and non-unitary cases. It is already known for a 
unitary case [57', that even though the amplitude of the triplet coupling is somewhat weakened with respect to the 
case of vanishing spin-orbit term, it is still finite. 

In Fig. [3] we show the anomalous Hall conductance in the superconducting phase for the two choices of the triplet 
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Figure 3: Anomalous Hall conductance for a spin triplet superconductor. Left panels present Hall conductance as a function 
of hz and \so for d = 1, whereas right panels as a function of d and X^o for hz = 0.5. Top figures correspond to the unitary 
case, see Eq. |11| while bottom figures correspond to the nonunitary case, Eq. |12[ 



pairing. We consider a unitary choice given by 

A^,-t- = d(-sinfcj^ +isinA;a;), A^^^ = A^_^ = 0, A^.^ = (i(sinfcj^ +?sinfc^) 

qx - 0, = 0, = 0. (11) 

and a nonunitary choice given by 

Af,^ = dsinfc^, A^^; = A^^^ = 0, A;^^ = 

qx = 0,gy = 0,gz = ySin^/Ca;. (12) 

In the case of unitary couphng (top panels of Fig. |3]), axy = if either Xgo — or hz ~ 0. However, in the case of a 
non-unitary couphng (bottom panels of Fig. |3]), axy = if Ago = 0, but for a nonzero spin-orbit coupling there is a 
finite Hall conductance even if hz = 0. In this nonunitary case there is a magnetization induced by the pairing, which 
leads to a finite a^y in a similar way as in '^He. 

In the unitary case the energy spectrum has a gap at the Fermi energy. This gap decreases as Xso increases. As Aso 
grows, the gap between the first and the second bands seems to decrease slightly and then it increases. In general, 
one can expects that small gaps between the bands will lead to large contributions to the Hall conductance. In the 
nonunitary case the energy spectrum also has a gap at the Fermi surface, which is small for small X^q, increases for 
slightly larger spin-orbit coupling, but vanishes when Xgo exceed Aso ~ 0.7. As Xso grows further the gap between the 
first and second bands increases. 

In the case when the superconductivity is intrinsic to the material, we have to solve the Bogoliubov - de Gennes 
equations self-consistently. We look for a situation of the type 

A^^-t- = J(-sinfc:r +sinfcj^), A^,; = A;_^ = 0, A;,; = (i(sinfca; -l-sinfcj^), 
qx = 0,qy = 0, = -^(-4sin/ca.sinfcy), (13) 



where the amplitude d is determined self-consistently for a given magnetization, taking into account that 

-y 



d= ^^(-sinA:^ + sinfcy)(V'g^V_fc^)- (14) 



where g is the pairing interaction. The corresponding numerical results are shown in Fig. |4j As the left panel 
shows, the superconductivity is destroyed for large enough spin-orbit coupling. In the right panel we see that the Hall 
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Figure 4: Spin triplet superconductor calculated self-consistently. Left: d as a function of and Xso- Right: Hall conductance 
along cuts of constant magnetization as a function of the spin-orbit coupling. 

conductance (as a function of A^o) decreases with increasing Aso, and as the transition to the normal phase appears, 
there are oscillations of the Hall conductance with relatively large amplitudes. 



III. AHE IN A CONVENTIONAL SUPERCONDUCTOR WITH MAGNETIC IMPURITY 

Consider now a classical spin immersed in a two-dimensional s-wave conventional superconductor. We use now a 
description of the system in the real space. In the center of the system, r = lc = {xc, Vc), we place a classical spin along 
the z direction. The kinetic energy part is described by a tight-binding model with hopping amplitude t, similarly as 
in the case of triplet superconductivity. The superconductor pairing is taken as s-wave, and the spin-orbit interaction 
[55] is assumed as in the preceding section. The electron operator is written in terms of the Bogoliubov operators, 

''Pir, c^) = X! ("«(^' - CTw„(r, ay-fl) . (15) 

n 

The zero momentum charge current in the ^ — x,y direction can be written as = l^f'^'f^^'^r^ where ■ijjf = 
( ipr,^- ^r,i )"^, and the velocity matrix operators are given by 

= J^{ltriU - KoTllcT^) . (16) 

Here f{r)f]1g{r) f(r + 5^)g{r) + f{r)g{r + 4) and f{r)-qtg{r) = f{r + S^)g{r) - f{r)g{r + 4)' where 4 is a 
displacement between nearest-neighbors along the /i direction, while ax,(Jy are Pauli matrices, as above. 
The real space wave functions obey the Bogoliubov - de Gennes equations for the energy excitations e„ 

/ -h-ep- JSpj-^ Xso{-Vx + iVy) \ 

At h + ep ~ JSpf^ XsoiVx - iVy) 

where h — tsg with sgf{f) = f{f+ 6). Furthermore, rjx = ±1 if the neighbor along a; is i^, 4- 1 (ix — 1) and r]y = ±1 
if the neighbor along y is iy + 1 (iy — 1). The parameter J describes coupling between the impurity spin and the 
spin density of conduction electrons. Note that the solution of this problem requires diagonalization of a AN x AN 
matrix, where N is the number of lattice sites. This is in contrast to the problem of triplet superconductor described 
in the previous section, where a partial diagonalization was possible due to the translational invariance. Owing to 
this symmetry, the problem could be reduced to a simple diagonalization of a 4 x 4 matrix for each momentum value. 
Since the effect of the magnetic impurity is rather local, a system of 15 x 15 lattice sites is sufficient to have small 
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Figure 5: In the top panels from left to right total magnetization, order parameter at the impurity site, and anomalous Hall 
conductance for the conventional superconductor with a magnetic impurity, calculated for a finite system including 15 x 15 
lattice points as a function of J and Xso- In the lower panel we show some cuts of the Hall conductance for J fixed and varying 
spin-orbit coupling. 



finite size effects, as we have shown previously [47J. We solve the problem self-consistently, as in a previous study (see 
Ref. Hi] for details). 

As in the case of triplet superconductivity studied in section II, the Hall conductance can be obtained from a Kubo 
like formula, which now reads 

- VS^.^^pVl.^^^-i5ur,{ri,a)*Vr,{r2,l)urn{ri,p)vra{r2,5)*^ (18) 

In this expression V is the complex conjugate, and Vp, means that the operator acts on the coordinate r^. 

The corresponding numerical results are presented in Fig. [5] where we show the total magnetization, order parameter 
at the impurity site, and the anomalous Hall conductance as a function of J and \so for a system of 15 x 15 lattice sites. 
The case with no spin-orbit coupling (A^o = 0) was studied before |1S]. Note, that the spin-orbit interaction shifts 
the critical value J^, at which the quantum phase transition to a magnetic state occurs, to higher values. However, 
if Xso is large enough, the transition is washed out. Note also that the various quantities reveal the quantum phase 
transition when fixing Xso and plotting them as a function of J. At the transition point, the impurity "captures" one 
electron and breaks a Cooper pair. Note that there is still a transition when we introduce the spin-orbit coupling, but 
one needs a larger coupling parameter J as the spin-orbit increases. If we increase the spin-orbit coupling further, 
superconductivity is destroyed and the Hall conductance exhibits strong oscillations, as in the case of the triplet 
superconductor . 

In order to emphasize the connection between behavior of the Hall conductance and the quantum phase transition, 
we show in Fig. |6] (for different spin-orbit couplings) the Hall conductance, amplitude of the order parameter at the 
impurity location, and the total magnetization of the conduction electrons as a function of the coupling between the 
spin density of the conduction electrons and the impurity spin. At the quantum phase transition both the amplitude 
of the order parameter and the total magnetization have discontinuities. At this critical coupling the Hall conductance 
has a sharp minimum which therefore signals the phase transition. 
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Figure 6: Order parameter at the impurity site, total magnetization, and anomalous Hall conductance for the conventional 
superconductor as a function of coupling strength to the impurity spin, calculated for the spin-orbit coupling corresponding to 
Xso = 0.8, \so = 1.1, and Aso = 1-4, as indicated. The first two values cross the quantum phase transition and for the highest 
value the transition turns into a crossover. 

IV. SUMMARY 

We have analyzed the anomalous Hall effect in superconductors, considering only the intrinsic mechanism that 
results from the interplay of Rashba spin-orbit interaction and magnetization. In the normal phase, the effect appears 
when both spin-orbit term and magnetization are nonzero. In a conventional spin-singlet superconductor of s-wave 
symmetry, an extended magnetization destroys the superconductivity. As we have shown, to have a nonvanishing 
anomalous Hall conductance in the superconducting phase it is then sufficient to assume a single magnetic impurity 
in the presence of spin-orbit interaction. However, vanishing coupling between the conduction electrons and magnetic 
impurity, or vanishing spin-orbit coupling, lead to zero Hall conductance. 

The case of a spin triplet superconductor is qualitatively different. An extended magnetization does not destroy the 
superconducting order. The magnetization, generally, either can be induced by an adjacent ferromagnet owing to the 
proximity effect (we may also consider the superconducting order as a proximity effect in heterostructures where some 
metal is coupled to a magnet and a superconductor) , or it may be an intrinsic property of the material (described by 
a self-consistent solution for the pairing amplitude). In the first case, two pairing forms lead to different results. If 
the pairing is unitary, the results are similar to those for the normal phase, and both magnetization and spin-orbit 
coupling are required for a finite Hall conductance. The superconducting case is very similar to the normal phase also 
when d is parallel to s. In the nonunitary case, however, there is a polarization associated with the pairing amplitude, 
and the Hall conductance is finite as long as the spin-orbit coupling is finite (a nonunitary pairing leads to a finite 
magnetization, as in the case of '^He). 

Since the spin-orbit interaction generates spin flips, its moderate values destroy superconductivity in both the 
conventional and triplet superconductors. In the case of s-wave superconductors, critical values of the spin-coupling 
J in the presence of spin-orbit coupling are larger than those for zero spin-orbit coupling, shifting thus the point at 
which the quantum phase transition appears. In the case of spin triplet superconductors with the pairing amplitude 
determined self-consistently, the spin-orbit coupling leads to suppression of the superconductivity through a continuous 
phase transition. Finally we have shown that the Hall conductance tracks the quantum phase transition induced by 
magnetic impurities in conventional superconductors. This provides transport measurement as a possible tool to 
detect the transition, related to earlier predictions that transport properties are affected by the presence of magnetic 
impurities in a superconductor |59]. We note that one of the interests of the AHE is that it can be easily measured. 

We acknowledge support from FCT grant PTDC/FIS/70843/2006. The work was partly supported by the Polish 
Ministry of Science and Higher Education as research projects in years 2011-2014 (VKD and JB). This work is also 
supported by the cooperation agreement between Poland and Portugal in the years 2009-2010. 
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